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Abstract: The Low Frequency Instrument (LFI) is a radiometer array covering the 30-70 GHz 
spectral range on-board the ESA Planck satellite, launched on May 14th, 2009 to observe the 
cosmic microwave background (CMB) with unprecedented precision. 

In this paper we describe the development and validation of a software model of the LFI pseudo- 
correlation receivers which enables to reproduce and predict all the main system parameters of 
interest as measured at each of the 44 LFI detectors. These include system total gain, noise tem- 
perature, band-pass response, non-linear response. The LFI Advanced RF Model (LARFM) has 
been constructed by using commercial software tools and data of each radiometer component as 
measured at single unit level. 

The LARFM has been successfully used to reproduce the LFI behavior observed during the LFI 
ground-test campaign. The model is an essential element in the database of LFI data processing 
center and will be available for any detailed study of radiometer behaviour during the survey. 

Keywords: Instruments for CMB observations; Space instrumentation; Microwave radiometers; 
Modeling of microwave systems. 
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1. Introduction 

The properties of the Cosmic Microwave Background (CMB) contain a wealth of information about 
physical conditions in the early universe and a great deal of effort has gone into measuring those 
properties since its discovery. Numerous ground based and balloon borne experiments, as well as 
space missions have been devised to obtain measurements of CMB spectrum, spatial anisotropics, 
and polarisation over a wide range of wavelengths and angular scales. Following the Cosmic Back- 
ground Explorer (COBE 1 ) and the Wilkinson Microwave Anisotropy Probe (WMAP 2 ) satellites, 
launched by NASA in 1989 and 2001 respectively, Planck 3 is an ESA space mission designed to 



http://lambda.gsfc.nasa.gov/product/cobe/ 
2 http://map. gsfc.nasa.gov/ 
3 http://www.esa.int/planck 



extract essentially all information encoded in the temperature anisotropics and to measure CMB 
polarisation to high accuracy. 

Planck, successfully launched on 2009 May, the 14th, will observe the entire sky from a Lis- 
sajous orbit around the second Lagrange point L2 of the Earth-Sun system. The focal plane formed 
by the 1.5 m off-axis telescope hosts two cryogenic instruments: the Low Frequency Instrument 
(LFI), covering 30-70 GHz in three bands Bersanelli et al. [2009], and the High Frequency Instru- 
ment (HFI) which covers the range 100 GHz to 857 GHz in six bands Lamarre et al. [2009]. The 
two instruments operate at 20K and at O.IK respectively. Such temperatures are achieved through 
a combination of passive radiative cooling and three active coolers Tauber et al. [2009]. While LFI 
and HFI alone have unprecedented capabilities, it is the combination of data from the two instru- 
ments that gives Planck the imaging power, the redundancy and the control of systematic effects 
and foreground emissions needed to achieve the scientific goals of the mission. 

The LFI has been designed to cover the low frequency portion of the Planck spectral range 
with three bands centred at 30, 44 and 70 GHz using pseudo-correlation receivers based on Indium 
Phosphide cryogenic high electron mobility transistors (HEMT) low noise amplifiers (LNA). 

In order to support the radiometers design and testing phase as well as the study the instru- 
ment behaviour and potential systematic effects during the survey we have constructed a software 
model of the receivers, the LFI Advanced RF Model (LARFM). In this paper we describe the devel- 
opment and validation of the LARFM, which has been constructed by using commercial software 
tools. The model was developed in two stages: the first stage, the Analytical LARFM (section is 
an analytical, parametric model built upon the requirements and specifications of the LFI receivers 
frequency by frequency. Its development started before the hardware was built and proved its use- 
fulness in predicting radiometers performance and to provide first-order estimates of the expected 
impact of potential systematics. 

The second stage is the Real-Data LARFM (section [l]), implemented after the LFI was built, 
which includes the data of each radiometer component as measured at single unit level channel by 
channel. The model is capable of reproducing in detail the main system parameters of each of the 
44 LFI channels "as built ", including total gain, noise temperature, bandpasses, non-linearity of 
the response. The Real-Data LARFM has been successfully used to reproduce the LFI behaviour 
observed during the LFI ground-test campaign Mennella et al. [2009a], Villa et al. [2009]. Cur- 
rently, the Real-Data LARFM is an essential element in the database of LFI data processing centre 
and is available for detailed studies of radiometer behaviour during the Planck survey. 

2. The Planck LFI design 

2.1 The overall LFI configuration 

LFI (Figure [I]) is an array of 1 1 horns, each feeding 2 orthogonal linearly polarised channels centred 
at 30, 44 and 70 GHz, on the focal plane of an aplanatic gregorian telescope. The LFI radiometers 
include the Front End Unit (FEU) and the Back End Unit (BEU), connected via 44 rectangular 
waveguides. Such a division into FEU and BEU has been performed in order to minimise the 
effects of power dissipation in the critical focal plane area. 

The Front End is cooled down to 20K by one of two Sorption Coolers, while the Back End, 
which is located in the body of the satellite, is at ambient temperature. The Front End Unit is the 
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Figure 1. A drawing of the LFI with its spacecraft interfaces. 



heart of the LFI instrument and it contains the feed array and associated OrthoMode Transducers 
(OMTs), hybrid couplers and amplifiers blocks. The Back End Unit comprises the Back End 
Modules (BEM) and the Data Acquisition Electronics. 

2.2 The LFI radiometer design 

The LFI array is made up of eleven Radiometric Chain Assembly (RCA) (Figure ||) each incor- 
porating the two radiometers connected to the feed horn. Each radiometer includes a pair of am- 
plification/detection chains correlated through a pair of hybrid couplers, constituting a continuous- 
comparison device Bersanelli et al. [2009]. In this scheme, the difference between the inputs to each 
of the chains (the signal from the telescope and that from a reference blackbody load, respectively) 
is continuously being taken. This scheme strongly suppresses 1/f noise generated by instabilities in 
the RF amplifiers, but not in the receiver elements following the second hybrid coupler. To remove 
this Back End 1/f, it is necessary to modulate the sign of the comparison using solid-state phase- 
shifters within the correlation section. The differencing receiver improves the stability at a cost of 
a factor of V2 in sensitivity compared with a total power scheme. Any difference between the sky 
and reference input signals can give rise to 1/f noise, in the LFI receivers this is further reduced by 
scaling one of the output data streams by a gain modulation factor "r". The blackbody reference 
itself must remain at a very stable temperature (see Valenziano et al. [2009]). 

The intensity AT of the smallest detectable signal is given by the radiometer equation (from 
Skou [1989]): 

Ar = ^ ZktZ^E , ( 2.i) 

VSeff • T 

where T m is the input temperature to the radiometer, r no i se is the radiometer noise temperature 
and t is the integration time. Radiometers are characterised by a spectral response g(v) through the 
radiometer band. The effective bandwidth is defined as: 
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Figure 2. The Low Frequency Instrument Radiometric Chain Assembly (LFI RCA) with the Data Acquisi- 
tion Electronics diagram. 
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3. The analytical LFI Advanced RF model 

The need to analyse the effect of LFI non-ideal response and to estimate the impact of systematic 
effects lead to the development of a radiometer system simulator Battaglia [2002], using the Agi- 
lent Advanced Design System (ADS) software. ADS provides Radio Frequency (RF) models for 
development of system specifications. It contains an RF simulator that predicts performance of 
complete RF systems and includes a set of block-level RF models for linear and non-linear compo- 
nents. The ADS Design Environment, Data Display, RF System Simulator and Models tools enable 
the graphic implementation of the receiver components and the possibility to obtain information 
about the RCA main properties, such as system gain and noise temperature. 
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3.1 Implementation 

The radiometer model is built as a network whose electromagnetic properties can be expressed 
using scattering parameters, taking into account the losses and reflections of the input signal, when 
it passes through each network element. 

For a generic multi-port network definition, it is assumed that each of the ports is allocated 
an integer n ranging from 1 to N, where N is the total number of ports. The insertion loss is the 
loss in load power due to the insertion of a component at some point in a transmission system. In 
measured insertion loss tables found in literature, the IS 21 1 2 is given rather than the pure insertion 
loss value. In order to evaluate the pure insertion loss, the measured return loss has to be taken into 
account. The IS21I is measured experimentally, in decibel, as 

\Sn\ 2 = ^. (3.1) 

'in 

52i includes both the effect of the insertion and the return loss, in fact: 



'IdB. I RL dB 

S 2 i = 10— Vl-10-nr, (3.2) 
and the pure insertion loss is given by: 

( IS21I 2 ^ 



IL dB = 10 Log w 
Return loss, instead, is given by: 



vi - 10 10 / 



(3.3) 



RL dB = \0Log w \S n \ 2 , (3.4) 
and total system gain is obtained by: 

G {dB) = 10 Log w ( — 1 , (3.5) 

where T m is the input load temperature and P out the corresponding power at detector output; Ay 
is the effective bandwidth and Kq is the Boltzmann's constant. r no i se , the system noise temperature 
is the input load temperature necessary to generate, in a noiseless system, the same output voltage 
of the noisy device with zero input load. 

3.1.1 Devices implementation 

Figure || shows a schematic of the implementation of the Analytic LARFM. The elements that 
constitute the RCA receiver shown in Figure ^ have been modelled in ADS starting from n-port 
device models, amplifiers, hybrids, wave-guides models and equation-based models. 
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Figure 3. The Analytical LARFM graphical representation in the ADS design environment. Each compo- 
nent is a port device completely characterised by the parameters displayed in the figure. 
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Input Loads. The receiver input load is generated by the port component termination. It is 
a 50Q resistance at temperature T ant that produces a power noise equal to P = KBT ant Av, where, 
in order to simulate the power absorbed by the sky and reference feed horns, T ant is the antenna 
temperature related to the thermodynamic temperature T by: 

hv 

Tant = -TT—T. (3.6) 
ew - 1 

Feed Horn and OMT. The sky and reference feed horns, as well as the OMT, are modelled 
as two-port devices, completely characterised by their scattering parameters - expressed in terms 
of insertion loss and return loss - and by their physical temperature, taking into account the noise 
temperature T no i se generated by the component itself into the network. Given an input signal T m at 
the input port of, e.g., a feed horn, the signal at the output of the element is given by: 

T ou t = \S 21 1 2 Ti n + T no i se . (3.7) 

Hybrids. The hybrid is modelled be a four port S -parameter set of equations using two values 
of magnitude and phase. They are set to correctly reproduce the in-phase and anti-phase addition 
of the inputs. The hybrid model by default has a perfect port isolation, but we tested also the effect 
of a crosstalk between the two paths of the signal (that simulates a non-ideal isolation between the 
two channels of the radiometer). The hybrid, being a noisy component, is also characterised by the 
physical temperature of the device. 

Front End and Back End Low Noise Amplifiers. The Low Noise Amplifier (LNA) is mod- 
elled using the amplifier component of the ADS object library. The LNA model uses scattering 
parameters for modelling the amplifier gain (5 21) and reflections at the input and output ports of 
the device (Sn and S 22); minimum noise figure (NF m i„) at optimum source reflection (S opt ) and 
equivalent noise resistance R n are used instead of noise figure to characterise the amplifier noise, 
providing also a correct return loss modelling response 4 . 

Phase Switch. The phase switch model is a two ports system. Scattering parameters are 
defined in polar coordinates, so that the change in phase is obtained by setting the S 21 parameter 
phase to 180 degrees. The temperature parameter characterises the noise properties of the phase 
switch model, as for the feed horn and OMT models. 

Waveguides. The RCA waveguides are modelled into two parts, the copper waveguide model 
and the gold-plated/stainless-steel waveguide model. Given the high thermal conductivity of cop- 
per, the copper waveguide model consists in a single rectangular waveguide element at the same 
temperature of the Front End Module, i.e. 22K. At the input and output ports two impedance trans- 
formers play a double role: they enable the waveguide matching to the characteristic impedance of 
the network (50 O) and they model the waveguide flanges return loss. 

The gold-plated and stainless-steel waveguide model, although characterised by the same param- 
eters, is more complex than the copper one, including several rectangular waveguide elements, 
taking into account the linear temperature gradient between the two flanges. 

In fact the stainless steel section interfaces to two heat sinks at either end, one at 300 K and the 
other at 22 K. The thermal distribution along the stainless steel waveguide is critical due to its high 

4 http://edocs. soco.agilent.com/display/ads2009/ Amplifier2+(RF+System+ Amplifier) 
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Figure 4. Schematic of the temperature distribution along the waveguides. V-grooves are thermal shields 
that isolate the cold Front End (22K) from the warm Back End Unit (293K). 



electrical resistivity. So, with reference to Figure f|, the stainless steel section has been modelled 
so that the temperature at the flanges and shields (V-groove) interfaces is preserved. For the gold 
plated sections the value for gold resistivity has been used. 

By setting the dimensions (width and height) of the waveguides section, the temperatures at 
the interfaces and the return loss at the flanges, the model is able to reproduce a cryogenic response 
of the waveguides, in terms of the insertion loss in the simulated frequency band. 

Bandpass filter. Its design is a Chebyshev response bandpass filter and its parameters (see 
Figure || for a complete list) completely constrain the insertion loss in the passband and the slope 
at the cut-off frequency. From the noise temperature analysis point of view, the band pass filter 
behaves like an attenuator, such as OMT and feed horns. 

Detector. The detector has been modelled with a two ports system characterised by a return 
loss value. The power conversion and integration over the band are performed by an appropriate 
calculation envelope, available in ADS, so that the output voltage integrated over the frequency 
band Av of the radiometer is given by: 

V out = f G-P out (v)dv, (3.8) 

where G is a gain factor and P ut{v) - V ou t(v) 2 /R, being V out (y) the simulated output voltage 
as a function of the frequency and R = 50 Q. the impedance of the network. 
Running ADS simulations it is possible to calculate: 

• scattering parameters over the band for the whole RCA model and component by component; 

• output voltage over the band; 

• integrated output voltage. 

3.2 Model verification 

The model is verified first at single component level and then at RCA level. 
3.2.1 Component level verification 

The functional consistency of each radiometer component is verified evaluating scattering parame- 
ters, output voltage and noise temperature. At single component level the S -parameters are verified 
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Figure 5. 30 GHz waveguide S 21 at ambient temperature (measured and simulated) and cryogenic tem- 
perature (simulated only). Ripples in measurements are due to reflections within the test set-up; those in 
simulations are due to the impedance mismatch at the flanges interfaces, since -30 dB return loss is set in 
this simulation. 



by the ADS primitive models themselves, so that a verification of sub-groups of cascaded elements 
has been carried out, e.g. cascaded feed horn and OMT. The simulated scattering parameters ex- 
actly match with those calculated using the scattering matrix theory, giving confidence that ADS 
applies the same expressions to compute scattering parameters for a cascaded system. 

Given the large quantity of cascaded elements, the validation of the waveguides model is 
performed by comparing the simulation results with the measured S 21 at room temperature. Then 
a prediction of their behaviour at cryogenic temperature is performed using the model, since no 
cryogenic measurements are available for the waveguides. Figure |5] show a comparison between 
simulation and measurements for the 30 GHz assembled waveguide (copper section and gold- 
plated/stainless steel section) at ambient temperature (300K). 

Cryogenic simulation shows that the S 21 is about ldB better than that measured and simu- 
lated at ambient temperature. The cryogenic model takes into account for the copper and gold- 
plated/stainless steel resistivity dependence on the temperature distribution along the waveguides. 



To evaluate the reliability of the model, functional tests and performance tests have been per- 
formed. Functional test goals included calibration of the system and a check of the correct mecha- 
nism for phase switching. 

The calibration factor, G, giving the conversion of the diode's output signal (in volts) into physical 
units (in kelvins), is defined as: 

V1-V2 AV 

where AT is the variation between two different input temperatures and AV is the correspond- 



ing variation in the output voltages. From equation it follows that the calibration factor G is 
independent of the temperature of the reference load. 

The following simulation set-up has been used to calculate the calibration factor: 
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Figure 6. Calibration test results at 30 GHz. The factor G is plotted in K/V versus frequency for fixed sky 
(left) and reference (right) load temperatures. The calibration factor is independent of phase switch status, 
labelled as {Lagi, Lagi). 



• Reference load temperature is fixed to 4K; 

• Sky load temperature is changed between 20K and 77K; 



All Phase Switches conditions have been analysed: both radiometer branches with degree 
lag; both branches with 180 degree lag; the intermediate status in which only one branch signal 
is 180 degree shifted (labelled (0, 180) and (180 , 0) in Figure §). Note that the effect of the lag 
on a single branch of the radiometer is to reverse the sky and reference signal at the outputs of the 
radiometer. The simulation has been repeated by fixing the sky load temperature to 4K and varying 
the reference load between 20K and 77K. Figure |6] shows the calibration factor G in the 30 GHz 
RCA frequency band, in any status of the two phase switches. Minor differences between the two 
results are due to asymmetries in the RCA (presence of the OMT in the sky input). 

3.2.2 System level verification 

Performance tests have been conducted, including system gain and system noise temperature cal- 
culus in the LARFM. All components parameters are set to their nominal value as from the Planck 
LFI specifications and requirements. The case of an RCA characterised by the insertion loss only 
is simulated; then simulations are repeated by taking into account for the reflections of the input 
signal (return loss). (Figure |7|) shows a generally good agreement with expected behaviour. Re- 
flections cause a general fall of the system gain value and the formation of band ripples, which 
contribute to an effective bandwidth reduction. System noise temperature is also affected by the 
presence of reflections at the waveguides flanges, causing ripples; moreover its mean value in the 
frequency band increases to about 1.5%. 

3.3 Applications 

3.3.1 Instrument performances with two FEM prototypes 

As a first application of the LARFM to the analysis of the instrument performances Franceschet 
[2004], we evaluated the system response in the frequency band in two different configurations 
including two prototypes of the Front End amplifiers at 44 GHz, whose performances, in terms of 
input/output return loss and total gain are measured in the frequency range 38-50 GHz (Figure [8]). 
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Figure 7. System Total Gain and System Noise Temperature of a 30 GHz RCA characterised by insertion 
loss only and both insertion loss and return loss. When reflections are calculated by the simulator, the total 
gain is reduced and the system noise temperature increases. 
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Figure 8. (Left) The first 44 GHz LNA prototype performances: input return loss, output return loss and 
gain measured at 293K. (Right) The second 44 GHz LNA prototype measurements at 295K. 



The LNA noise figure is given by the measurements on a third LNA and has not changed during 
the entire analysis of the prototypes; its profile is visible in System Noise Temperature simulation 
for the nominal case (see Figure |9[ right). 

In order to evaluate the effect of the Front End LNA performance on the system total gain, 
noise temperature and effective bandwidth, we set all RCA parameters, with the exception of those 
of the Front End amplifiers, to their nominal value at 44 GHz. The inland performances of the 
two LNA prototypes are measured at ambient temperature (a worst case). The nominal case, with 
the LNA input return loss RL in = -5 dB, output return loss RL uut = -6 dB and gain G = 30 dB, is 
also simulated for a comparison. 

Figure |9] shows that in both simulations including LNA prototypes, the total gain follows the 
shape of the corresponding prototype, with wide ripples for frequencies above 47 GHz, due to the 
effect of the input return loss. As expected, the RCA with the second LNA prototype exhibits better 
performance in terms of total system gain, in particular for frequencies above 44 GHz. Simulations 
show similar results for the system noise temperature: it rapidly diverges for frequencies above 48 
GHz in the first prototype simulation; the input return loss of the LNA plays a leading role in this 
simulation, since nearly all incoming power is reflected at the LNA input port. 

The effective bandwidth, as calculated by the LARFM for all test cases is summarized in 
Table |. 
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Figure 9. System Total Gain and Noise Temperature simulated for 44 GHz RCA's including two LNA 
prototypes and the nominal case (LNA performances as from requirements). 



Frequency Range 


38-50 GHz 


39.6^8.4 GHz 


B e ff Prototype 1 


6.93 


6.28 


B e ff Prototype 2 


8.94 


7.45 


B e ff Nominal 


11.75 


8.84 



Table 1. Effective bandwidth B e ff in the measured data frequency range (38-50 GHz) and in the Planck 
LFI required frequency bandwidth (20% of the centre frequency). 



This analysis shows that the LARFM is a useful tool to investigate the behaviour of the Planck 
LFI instrument when in band non idealities characterise the system components. 

3.3.2 Effects of OMT asymmetry 

If the OMT arms have different performance, in terms of loss and reflection in the frequency band, 
the two orthogonal components of the polar signal are affected differently, and will be confused 
with a real polarisation signal when processing the polarisation information. As an application of 
the LARFM, the impact of asymmetries, in terms of losses and reflections, in each one of the OMT 
arms on the in-band total system gain is studied. Simulations are performed by setting all RCA 
parameters to their nominal values, with the exception of the OMT insertion loss and return loss, 
measured on the Qualification Model (QM) of the 30 GHz OMT. 

Simulation results in Figure [10| shows that the maximum gain mismatch between the arms of 
the radiometer is about 0.25% at 30 GHz. 

For a description of the impact of asymmetries between the two arms of the 30 GHz RCA on 
the CMB polarisation measurement, the following expression has to be evaluated: 

A (AT) = ±\f Gi (v)T sky (v)dv- J G l (v)nT ref (v)dv] + 

-£UG2(y)T sky (v)dv-jG 2 (v)r 2 T ref (v)dv] ' 

where G\{v) and G 2 {v) are the gain relevant to the two arms of the radiometer (Figure |To| ), 
T s ky and T re f are the observed sky temperature and the reference load temperature and r\ and r 2 
are the gain modulation factors for both arms of the RCA. The sky signal is parametrised assuming 
a power-low behaviour for the sky temperature, 
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Figure 10. Effect of the OMT in band asymmetries on the total system gain of the 30 GHz RCA Qualification 
Model. 
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Figure 11. Difference, in mK, between the final output of the two branches of the radiometer. 



T sk y = a V -' s (3.11) 

where a and are derived from WMAP foreground maps at 23, 33, 41, 61 and 94 GHz and 
the CMB signal expression, 

hv/K B T 

rcm(y) %w,-i To (3 - 12) 

where To = 2.725 K, and Kb is the Boltzmann's constant. The K, Ka and Q band data fix the a 
and /3 values in the ranges a € [5 • 10 3 ,6 • 10 3 ] and /? € [0.25,0.35]. The reference load temperature 
follows the frequency dependant 4.8K black-body curve for this analysis. 

By calculating the integrals for the specified variation ranges of the normalisation parameter a and 
of the spectral index /?, the difference between the two differenced output signals, i.e. the expression 



A(AT), is plotted in Figure |11 



4. The Real-Data LFI Advanced RF Model 

During the qualification and performance tests, each component of the LFI radiometers was inde- 
pendently characterised in terms of frequency response. These hardware measurements replaced 
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the parameters used in the analytical version of the LARFM with the objective of estimating the 
effective response of the 44 channels of the LFI flight hardware. 

The new LARFM maintains the same implementation just for the waveguides simulator, while 
the other units are replaced by their measurements. Moreover, the main objective of the LARFM 
became the estimation of the bandpass response, see Zonca et al. [2009], due to the high impact 
of systematic effects on end-to-end bandpass measurements. The next sections explain the imple- 
mentation and the results obtained with the LARFM. 



4.1 Implementation 



Figure [12J shows a schematic of the implementation of the Real-Data LARFM. In this model the 
pseudo-correlation strategy is not considered, see section |4.1.2| , so each channel is completely 
independent from the other channel of the same radiometer. Moreover, the number of devices is 
strongly reduced because FEM and BEM devices where fully characterized after assembly. 
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Figure 12. Model schematic for each LFI channel: feed horn and OMT are characterized by their S- 
parameters; FEM by gain and noise temperature; waveguides are simulated analytically; BEM by return 
loss, transfer function and integrated noise temperature. 



4.1.1 Feed horn and OMT 

Feed horn and OMT return losses were measured on their assembly, injecting the test signal from 
the OMT port to be connected to the FEM. Insertion losses, instead, are available for the OMT only, 
but the impact of the feed horn is negligible. Feed horn and OMT model is a single Sparameters 
passive component based on these data. 

4.1.2 FEM 

A FEM model considering the pseudo correlation strategy needs spectral response measurements 
of each component assembling the FEM: 

• Hybrids 

• Low Noise Amplifiers 

• Phase switches 

• Internal FEM waveguides 
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However, all these components when tightly assembled into the FEM change their response 
significantly due to mutual interactions. It is therefore more reliable to use the data from the tests 
of the complete FEM. 

Each channel is simulated independently from the others using data from gain and noise tem- 
perature tests performed by the FEM producers. Each channel is considered as a total power ra- 
diometer frozen in the configuration of the phase switches states which connects it to the sky signal. 
Indeed, there are two phase switch state combinations where the same BEM channel is connected 
to the sky arm, but, after phase switch balancing, the responses match better than ~ AdB. 

In ADS, the FEM has been modelled as an amplifier with gain given by the performance tests 
on each FEM channel and noise figure computed by the system temperature as: 

FEM S 11 could only be measured at ambient temperature, therefore the LARFM implements 
the mean expected value 5 n = -IdB. The RCA is then studied in an ideal configuration, in which 
all four RCA channels are looking simultaneously to the sky signal. 

4.1.3 BEM 

The BEM model is made up of: 



• Reflection only Sparameter component for return loss 

• A noise generator for the noise temperature 

• An equation object for the gain 



In order to model the whole RCA, the RF signal coming from the FEM through the waveguides 
is summed with a noise generator based on the measured BEM noise temperature. This signal is 
then amplified by the measured spectral transfer function converting input power to output volts at 
a given frequency, which includes the effect of all the BEM components together. 

The output is a bandpass response of the complete RCA which can then be integrated over 
frequency in order to compute the detector voltage output that can be compared with test results. 

During the flight hardware test campaign, 30 and 44 GHz RCAs showed signal compression 
at BEM level, i.e. the BEM output doesn't increase linearly with the input temperature. This 
effect has a big impact on system temperature calculation, because the data extrapolation process 
emphasises even a small non-linearity. Therefore, a parametric non linear analytic model was 
developed in Mennella et al. [2009b] in order to fit test data for extracting the overall RCA gain, 
system temperature and BEM compression factor. The model is based on the following equation: 



V ut = Go ■ 



1 



■(T A + T noise ), (4.2) 



1 + b • Go • (T A + T noise) _ 

where Ta is the antenna temperature, Go is the linear gain and b is the compression factor. The 
compression factor is a single scalar value per channel which represents the BEM gain dependence 
on the input power. BEM tests were performed at a power level where compression is expected. 
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The first operation to perform then is to calculate the linear gain by using the known information 
on the compression factor measure in the relevant tests. The equation used is: 



Glin - 



G t , 



(4.3) 



1 - b ■ P in ■ G test ' 

where G tes t is the total BEM gain or transfer function (i.e. output voltage vs input power), b 
is the scalar compression factor, P\ n is the input power for the BEM test and Gu„ is the computed 
linear gain. The last step is to substitute Gu„ for the gain parameter of the BEM and change the 
equation describing the BEM to: 



GP h 



(l+b-P in -G) 



(4.4) 



4.2 Results 



The Real-Data LARFM outputs are the channel bandpasses, volt outputs, gains and noise temper- 
atures either on the full RCA or in any section of it. These data are the basis of several analyses on 
the LFI radiometers, in the following paragraphs we are going to introduce the most interesting. 
The most important result of the LARFM is bandpass estimation, see Zonca et al. [2009]. 



Figure |13| shows a good agreement between LARFM bandpasses and end-to-end measurements, 
the end-to-end measurements are however affected by systematic effects and their quality is not 
sufficient for a Quantitative analysis. 
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Figure 13. LFI27S-1 1 at 30 GHz. Comparison of LARFM bandpasses with end-to-end swept source mea- 
surements performed on WR28 nominal bandpass, between 26.5 and 40 GHz, see Zonca et al. [2009]. 



Figure [14] shows the output voltage as a function of input load temperature measured in a 
linearity test compared to the LARFM results assuming a linear BEM response and implementing 
the BEM non-linear model as explained in section [4.1.3[ The compression parameter b of Eq. 4A 
was measured using the test results and the comparison shows a good agreement between the non- 
linear model and the measurements. 
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Figure 14. Output voltage versus temperature for a linearity test compared to the LARFM linear and non- 
linear BEM models. 



The LARFM was also successfully applied to predict the impact on bandpasses of the sub- 
stitution of a damaged FEM unit with a spare unit, to assess whether it would be acceptable to 



replace or necessary to repair the unit. Figure 15 shows the impact on the LARFM bandpasses of 
the eventual substitution of the flight FEM with the spare unit on RCA 24. The response of the 
spare FEM is more structured than the nominal unit and the bandwidth is reduced by about 20% 
(see Zonca [2009]), the unit was successfully repaired. 



FS vs FM SIMULATIONS (10 log mV) 




Figure 15. Comparison between LARFM bandpasses with nominal versus spare FEM on LFI24M-00 



The RF model has been used also during the cryogenic test campaign to successfully investi- 
gate the saturation at the BEM diode, by modelling the RF power level at the diode input. 

A more direct use of the model is foreseen in the unfortunate case of the radiometers working 
in non-nominal configurations, for example different thermal conditions. In this case it is possible 
to estimate the impact of these non-nominal conditions on each component and then make use of 
the LARFM to simulate the non-nominal behaviour in terms of bandpasses and output voltage. 
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5. Conclusions 



We have described the LFI Advanced RF Model, based on commercial software tools and measure- 
ments of the radiometer component at single unit level. The LARFM provides a good description 
of the response of each of the 44 channels of the Planck-LFI array, such as radiometer response 
(including non-linear behaviour), system gain, noise temperature, spectral shape of the radiometers 
bands. 

The measurements of the scattering parameters of the single units (i.e., feed horns, Ortho- 
Mode transducers, front-end modules, waveguides, back-end modules) are very precise (typically 
better than O.ldB), while in some cases the end-to-end measurements of the integrated radiometer 
properties are more difficult. This is the case, for example, for the evaluation of the shape of the 
spectral bands, whose knowledge is particularly important for polarisation analysis. In these cases, 
the availability of the LARFM has proved particularly useful. 

After being successfully used to support the Planck-LFI design and ground-test campaigns, 
the model is now incorporated in the LFI data processing centre to allow detailed studies of ra- 
diometer behaviour during the Planck survey. Similarly, future precision radiometric instruments 
will likely require similar models to be developed to support instrument development and data 
analysis. Our experience has shown the usefulness of such an RF software model and indicates a 
strategy and some technical tools that may be useful for other projects in the future. 
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